Odd triplet superconductivity in clean and moderately disordered SFFS junctions 
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We study the Josephson effect and pairing correlations in SFFS junctions that consist of con- 
ventional superconductors (S) connected through two metallic monodomain ferromagnets (F) with 
transparent and spin inactive interfaces. We solve the Eilenberger equations for arbitrary relative 
orientation of magnetizations of the two F layers in the clean limit and for moderate disorder in fer- 
romagnets. Spatial variation of pair amplitudes, singlet f s , and odd in frequency triplet fto and fti, 
with and ±1 spin projections, as well as the Josephson current-phase relations are calculated for 
different values of the ferromagnetic layers thickness and angle a between in-plane magnetizations. 
In contrast to the dirty limit case, we find that for < a < -k both spin singlet and triplet pair 
amplitudes in F layers power-law decay in the same oscillatory manner with distance from the FS 
interfaces. This decay gets faster as the impurity-scattering rate in ferromagnets is increased. The 
computed triplet amplitude fti has opposite signs in the two magnet regions, penetrates into the 
superconductors and monotonically decays over the same distance, which is the superconducting co- 
herence length, as the singlet amplitude f s saturates to the bulk value. We point out that influence 
of misorientation of magnetizations on the Josephson current can not be attributed directly to the 
appearance of odd triplet correlations. 

PACS numbers: 74.45.+C, 74.50.+r 



I. INTRODUCTION 



Odd-frequency triplet superconducting correlations 
with non-zero total spin projection induced by proxim- 
ity effect in heterostructures containing superconductors 
with ordinary singlet pairing and inhomogeneous ferro- 
magnets have been of considerable interest in the last 
decade Odd-frequency pairing mechanism was 
proposed long ago in an attempt to describe the A phase 
of superfluid 3 He. 7 Although this superconducting state 
is thermodynamically stable,— it was found that in the 
case of superfluid 3 He the pairing is odd in space (p- 
wave) rather than in time, "odd" triplet superconduc- 
tivity occurs in certain superconductor-ferromagnet (SF) 
structures.— It is believed that this exotic s-wave pairing 
state which is even in momentum, but with the triplet 
correlations being odd in frequency, can enhance signif- 
icantly the superconducting penetration length at the 
SF interface^ 3 - providing a Josephson current through 
half-metallic barriers4££ The simplest examples of SF 
heterostructures with inhomogeneous magnetization are 
FSF and SFFS heterojunctions with homogeneous mon- 
odomain ferromagnetic layers having an angle a between 
their in-plane magnetizations. These structures have 
been studied using quasiclassical approach in both dif- 
fusivei^iSiiSiiiiiHii 3 . and cleanii^ii^ limits by solving 
Usadeli 7 - and Eilenberger— equations, respectively, and 
by solving the Bogoliubov-de Gennes equation i 6 ' 19 i 20 i 21 i 22 
In the case of parallel (a = 0) magnetizations, SFS 
junctionS ) 23 i 24 besides spin-singlet superconducting cor- 
relations, only odd triplet correlations with zero total 
spin projection exist. These correlations penetrate into 
the ferromagnet over a short length scale determined 
by the exchange energy. For noncollinear magnetiza- 
tions, triplet correlations with nonzero spin projection 
are present as well. It is expected that they are not sup- 



pressed by the exchange interaction, and consequently 
that they are long ranged^ It has been predicted for 
diffusive junctions that long-range spin-triplet compo- 
nents with spin projection ±1 should have a dramatic 
impact on transport properties and the Josephson effect, 
displayed through a nonmonotonic dependence on an- 
gle between magnetizations. In diffusive Josephson junc- 
tions, the length scales associated with short- and long- 
range correlations are, respectively, £ F = ^JhDp/ho and 
£jy = \JhDp /ksT, where Dp is the diffusion constant 
in the ferromagnet, and thermal energy kpjT is typically 
much smaller than the exchange energy ho. However, this 
is not the case in ballistic SF heterostructures in which 
the only characteristic length is the ferromagnet coher- 
ence length £p = Hvp/ho, where vp is the Fermi velocity. 
In the clean limit, influence of misorientation of magneti- 
zations on the Josephson current or conductance can not 
be attributed to the triplet correlations. 1 ^ 20 Moreover, 
for noncollinear magnetizations, a^O and 7r, both spin 
singlet and triplet pair amplitudes in ferromagnetic lay- 
ers decay in the same oscillatory manner with distance 
from the FS interfaces. 21 ' 22 

In this paper, we study the Josephson effect and odd- 
frequency triplet superconductivity in clean and moder- 
ately disordered SFFS junctions where the magnetic in- 
terlaycr consists of two monodomain ferromagnets hav- 
ing relative angle a between their in-plane magnetiza- 
tions. Solving the Eilenberger equations we calculate 
singlet and triplet pair amplitudes, and the Josephson 
current. In contrast to the dirty limit case we find that 
for < a < 7T, both spin singlet and triplet pair ampli- 
tudes in F layers are oscillating and power low decaying 
with distance from the FS interfaces. This decay gets 
faster as the impurity-scattering rate in ferromagnets is 
increased. The computed triplet amplitudes have oppo- 
site signs in the two magnet regions, penetrate into the 
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superconductors and monotonically decay over the su- 
perconducting coherence length, = hvp/irA , where 
Ao is the superconducting pair potential, which is the 
same distance on which the singlet amplitude saturates 
to the bulk value. As in the previous workf^ no substan- 
tial impact of spin-triplet superconducting correlation on 
the Josephson current has been found. The critical value 
of the Josephson current I c is a monotonic function of 
angle a when the junction is far enough from the — 7r 
transitions. Otherwise, J c (a) is a nonmonotonic function 
of a with characteristic dips related to the onset of — it 
transitions, but this nonmonotonicity is a consequence of 
decreasing influence of the exchange potential with in- 
creasing misorientation of magnetizations,^ 9 Influence of 
misorientation of magnetizations on the Josephson cur- 
rent can not be attributed to the appearance of triplet 
correlations: with increasing misorientation of magneti- 
zations, a —> 7r, the Josephson currents are the same as 
for parallel magnetizations, a — 0, with correspondingly 
decreasing exchange potential, ho —>■ 0. Therefore, ro- 
tation of magnetizations in SFFS junctions is equivalent 
to the decreasing influence of the homogeneous exchange 
potential, in the absence of triplet correlations. 

The paper is organized as follows. In Sec. II we 
present the model and the solutions of Eilenberger equa- 
tions that we used to calculate the Josephson current and 
spin-singlet and -triplet pair amplitudes. In Sec. Ill we 
present the numerical results, and the conclusion in Sec. 
IV. 



II. MODEL AND FORMALISM 
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FIG. 1: Schematics of an SF1F2S heterojunction. The mag- 
netization vectors lie in the y-z plane and form the opposite 
angles ±q/2 with respect to the z-axis. 



subara frequencies w n — TrkBT(2n + 1), n = 0, ±1,±2, 
etc. Spin indexes are a =|, J.. 

The Eilenberger equation in spin<g>particle-hole space 
can be written in the compact form 

hv x d x g + u> n f 3 ® 1 - iV + A + h(g)/2r,g = 0, (2.1) 

with normalization condition g 2 = 1. We indicate by 
• • • and • : • 2 x 2 and 4x4 matrices, respectively. The 
brackets (...) denote angular averaging over the Fermi 
surface (integration over 9) and [. . .] denotes a commu- 
tator. Here, the quasiclassical Green functions 
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(2.2) 



We consider a clean SF1F2S heterojunction consist- 
ing of two superconductors (S), and two uniform mono- 
domain ferromagnetic layers (Fi and F2) of thickness d\ 
and g?2, with misorientation angle a between their in- 
plane magnetizations (see Fig. [1]). Interfaces between 
layers are fully transparent and magnetically inactive. 
Superconductors are described in the framework of qua- 
siclassical theory of superconductivity, while the ferro- 
magnetism is modeled by the Stoner model, using an 
exchange energy shift 2ho between the spin subbands. 
Disorder is characterized by the mean free path I — vft, 
where r is the average time between scattering on im- 
purities and vf is the Fermi velocity assumed to be 
the same everywhere. We consider the clean limit and 
moderately diffusive case when I is larger than the two 
characteristic lengths: the ferromagnetic exchange length 
£f = hvp/ho, and the superconducting coherence length 
£s = Tiuf/VAo, where Ao is the superconducting pair 
potential. 

In this model the Eilenberger Green functions 
9aa'( x ^ v x,^n), f aa > (x, v x , u n ), and f+ a ,(x,v x ,uj n ) de- 
pend on the center-of-mass coordinate x along the junc- 
tion axis x, on the angle 9 of the quasiclassical trajec- 
tories with respect to x axis, v x — vpcos9 being the 
projection of the Fermi velocity vector, and on the Mat- 



are related to the corresponding Gor'kov-Nambu Green 
functions G integrated over the energy £k = h^k 2 /2m — /i, 



9 = ~ I de^G. 



(2.3) 



The matrix V is given by 



V 



Re h(x) 



if 3 ® Im h(x) ■ <x , (2.4) 



where components a x ,a y , a z of the vector <r and fi , T2 , T3 
are Pauli matrices in spin and particle-hole space, respec- 
tively. The in-plane, y-z, magnetizations of the neighbor- 
ing F layers are not collinear in general, and the magnetic 
domain structure is described by the angle Tct/2 with 
respect to the z-axis, for the left (Fi) and right(F 2 ) fer- 
romagnets, respectively. For simplicity, we assume equal 
magnitude of the exchange interaction in ferromagnetic 
domains, h(a;) = /io(0, T sin a/2, cos a/2), and A(x) = 
for — di < x < d2- 
The matrix A is 



A = 



fiA 
fiA* 



(2.5) 
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We assume that the superconductors are identical and, 
in the step-wise (non self-consistent) approximation, we 
take the pair potential A in the form 



A = An 



e -it/2 Q (_ x _ di j + e iH*®(x _ d2 )l ( (2.6) 



where Ao is the bulk superconducting gap and <j) 
is the macroscopic phase difference across the junc- 
tion. The temperature dependence of Ao is given by 

A (T) = A (0)tanh(l.74^T c /T- l). 25 Although the 

self-consistent calculations are needed when the proxim- 
ity effect is strong between S and F layers^ one would 
not expect this to qualitatively change the spatial varia- 
tion of the Green functions. 

In the clean limit, t — > oo, solutions of Eq. (|2.ip for 
the left superconductor (x < — g?i) can be written in the 
usual form for normal Green functions 



9u( x ) 




(2.7) 


9n( x ) 


= D 2 e^ x , 


(2.8) 


9ll( x ) 


= D 3 e K ° x , 


(2.9) 




- -FT + U ^ e 


(2.10) 



with 



2uj n + hv x K s 
A 2A 



and for anomalous Green functions 
/TtW ^ " 



fi„ 2uj n + HVxKi 

A 2A 



2A 

2u) n + hv x n. 



+ flV x K s 

■D 3 e K * x , 



D ie KsX , 
D 4 e K ° x , 



Hvf \ cos(#) 



(2.11) 
(2.12) 
(2.13) 
(2.14) 



(2.15) 



and n„ = y/ oj 2 + | A | 2 . For the right superconductor 
(x > d 2 ) 1 the solutions keep the same form with k s — ► 
— k s and new set of constants D 1 , D 4 . 

Solutions for the Green functions in the clean limit, 
t — > co, for the left ferromagnetic layer Fi, — di < x < 0, 
can be written in the form 



and 



9u( x ) 
911 0) 



K! + it&n(a/4)K 2 e iK0X - it&n{a/4)K 3 e- lK0X , 
K 2 e lK ° x + t&n 2 {a/4)K i e- lK0X + - t(m(a/2)[K 1 - K 4 



= K 3 e 



tan 2 (a/4)if 2 e i 



tan(a/2)[Xi - K 4 



K 4 -ita,n(a/4)K 2 e iKoX + i tan(a/4)fsT 3 e" 



f n (x) = C ie - KX -it&n(a/4)C 2 e- K - x ~ itnn(a/A)C 3 e~ K + x , 
f n (x) = C 2 e- K - X - tan 2 (a/A)C 3 e- K+x - i tan(o/2)[Ci + C A ]e~ K 

f L1 (x) = -C 3 e~ K + x + tSin 2 (a/A)C 2 e- K - x + % - tan(a/2)[Ci + C 4 ]e 
f u (x) = C A e~ KX -it&n{a/A)C 2 e- K - x - it&n(a/A)C 3 e' K + x , 

I 



(2.16) 
(2.17) 

(2.18) 
(2.19) 

(2.20) 
(2.21) 

(2.22) 
(2.23) 



where 



and 



k = 2h /hv x , 



K = 2uj n /hv x , 
K- = 2(uj n - ih )/hv x , 
k+ = 2(u n + iho)/hv x . 



(2.24) 



(2.25) 
(2.26) 
(2.27) 



The solution for the right ferromagnetic layer F 2 , < 
x < g?2, can be obtained by substitution a — ► —a, with 



a new set of constants Ki,...,K 4 and C 1 ,...,C 4 . The 
complete solution requires one to determine 24 unknown 
coefficients: 4+4 in superconducting electrodes, and 8+8 
in ferromagnetic layers. Boundary conditions of conti- 
nuity of g at the three interfaces provide necessary 24 
equations. 



For ho = solutions (|2.16[) - (|2.19j) reduce to the solu- 
tions for SNS junctions with K 2 = K 3 = and K 4 = K\. 
Solutions (j2~2T)l) - (f2~2"3"|) reduce to the SNS case if we take 
a = and than ho = with C% = C4 = and C 3 = C 2 . 
Note that the same solutions are obtained for a = it and 
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The supercurrent density is given by the normal Green 
function through the following expression 

7(0) = 7reiV(0)A; s T^^(t; x Im g^(v x )), (2.28) 

where N(0) is the density of states per spin at the Fermi 
level. 

The zero-time pair amplitudes, singlet f s , and triplet 
/to and fti, with and ±1 projections of the total spin of 
a pair, are defined in terms of anomalous Green functions 
as 

/, = nN(0)k B Tj2(fn ~ /it), ( 2 - 2 9) 

f t0 = nN(0)k B Tj2(fn + fnl (2-30) 

fti = nN(0)k B Tj2(f n + f lt ). (2.31) 

The ground state in SFS junctions can be or ir stated 
For <fi — the singlet amplitude f s is real, while the 
triplet amplitudes f t Q and f t \ are imaginary, and for 4> = 
7r the opposite is true. For < <f> < 7r, all the amplitudes 
are complex functions. In the figures, we will normalize 
all amplitudes to the value of f s in bulk superconductors, 

f 8b =7rN(0)k B Tj2^- (2-32) 

The influence of moderate disorder in ferromagnetic 
layers on the Josephson current and on the pair am- 
plitudes is studied in the limit hor/h 3> 1. From the 
term h(g)/2r in Eq. (12. ip we kept the largest contribu- 
tion, which is (<7<t<t) calculated for r — > oo. The contri- 
butions of {f aa ') and (<7ctct) are neglected. Introducing 
(.9Tt) = R + H, we find that (gn) = R — il. Numerical 
calculations show that R is independent of x, and R ks 
sign(w n ) for any a. The imaginary part I oscillates with 
x, except for a = 0, but it is nontrivial only in the vicin- 
ity of to n = 0. Therefore, the contribution of I to the 
pair amplitudes, as well as to the supercurrent, can be 
neglected, and the previous solutions Eqs. (|2.16p - (|2.23p 
are used by replacing only k, k±, Eqs. (|2.25|) — (|2.27|1 . with 
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FIG. 2: (Color online) The current-phase relation I((f)) in 
the clean limit for T/T c = 0.1, h /E F = 0.1, di = d 2 = 
50/cp 1 , and five values of the misorientation angle: a — 
(thick solid curve), n/3 (dotted curve), n/2 (thin solid curve), 
2n/3 (dashed curve) and n (dash-dotted curve). The latter 
coincides with I{4>) in the corresponding SNS junction (ho — 
0). 
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FIG. 3: (Color online) The current-phase relation I((f)) in 
the clean limit for T/T c = 0.1, h /E F = 0.1, di = d 2 = 
500/cp 1 , and five values of the misorientation angle: a — 
(thick solid curve), 7r/3 (dotted curve), 7r/2 (thin solid curve), 
27r/3 (dashed curve) and n (dash-dotted curve). The latter 
coincides with I(<j>) in the corresponding SNS junction (ho = 
0). 



k = kH , (2.33) 

R 

k_ = k_ H , (2.34) 

k+ = k.+ + . (2.35) 

III. RESULTS 

For simplicity, we consider a single transfer channel 
(one-dimensional case) . We illustrate our results for sym- 
metric junctions with relatively thin and thick ferromag- 
netic layers, di = = 50k p 1 and 500/c^ 1 , and for low 



temperature T/T c — 0.1. Superconductors are char- 
acterized with the zero temperature value of the bulk 
pair potential Aq(0)/-E'f = 10~ 3 , which corresponds to 
£s(0) = 636A;^ 1 . In order to achieve quantitatively reli- 
able results within the quasiclassical approximation, we 
assume that all interfaces are fully transparent, the Fermi 
wave vectors in all metals are equal (&f), and the ferro- 
magnets are relatively weak, ho/Ep = 0.1, which corre- 
sponds to £p = 20fc^ 1 . 

Ballistic regime, t — ► oo, is illustrated in Figs. [2][5] 
The Josephson current is calculated from Eq. (|2.28[) for 
a single transverse channel (the angular averaging (...) 
reduces to one half of a discrete sum over 9 — and 
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FIG. 4: (Color online) Spatial dependence of singlet and 
triplet pair amplitudes f a , fto and fti, normalized to the bulk 
singlet pair amplitude f a b, for T/T c = 0.1, Ho/Ef = 0.1, 
d\ — di — 50/cp 1 , and three values of the misorientation an- 
gle: a — (thick solid curve), tt/2 (thin solid curve) and tt 
(dash-dotted curve). All amplitudes are calculated for (f> — 
and = 0. 



FIG. 5: (Color online) Spatial dependence of singlet and 
triplet pair amplitudes f a , fto and fti, normalized to the bulk 
singlet pair amplitude f s b, for T/T c — 0.1, Ho/Ef = 0.1, 
di — d,2 = 500/cJ 1 , and three values of the misorientation an- 
gle: a — (thick solid curve), tt/2 (thin solid curve) and tt 
(dash-dotted curve). 



8 = tt, and one-dimensional density of states per spin 
is iV(0) — 1/2ttHvf)- We illustrate the current-phase 
relations I(cf>) in Figs. [5] and [3] for five values of the 
angle between magnetizations in ferromagnetic layers, 
a = 0, 7r/3, tt/2, 27r/3, and tt. Previous results obtained 
in Ref. 19 by solving the Bogoliubov-de Gennes equation 
are exactly reproduced for this case of transparent inter- 
faces and equal Fermi wave vectors. In particular, we 
see in Fig. [2] that the transition between and tt states 
can be induced by varying the misorientation angle a, in- 
cluding the coexistence of stable and metastable and tt 
states with dominant second harmonic. 27 However, non- 
monotonic dependence of the Josephson current I{4>) on 
a is due to decreasing influence of the exchange potential 
with increasing misorientation of magnetizations: I(4>) 
for noncollinear magnetizations (a ^ 0) is practically 
equal to the one for homogeneous magnetization (a = 0) 
with correspondingly smaller exchange energy. Influence 
of a on /(</>) can not be attributed directly to the appear- 
ance of odd triplet correlations. When the junction is far 
enough from the — tt transitions, the critical Josephson 
current is a monotonic function of a, Fig. [3] We em- 
phasize that I(4>) curves for a = tt are the same as in 
the corresponding SNS junctions (ho = 0) since the in- 
fluence of opposite magnetizations in F layers of a sym- 
metric SFFS junction cancels out in the clean limit and 
for transparent interfaces. 

The spatial variation of the pair amplitudes is shown 
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FIG. 6: (Color online) Spatial dependence of singlet and 
triplet pair amplitudes / s , fto and fti, normalized to the bulk 
singlet pair amplitude f s b, for T/T c — 0.1, ho/E^ = 0.1, 
di = d2 — 50/cp 1 , and two values of the misorientation an- 
gle: a — (thick solid curve) and tt/2 (thin solid curve). All 
amplitudes are calculated for <j> = tt and 6 = 0. 
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FIG. 7: (Color online) Real R and imaginary I parts of (<?n) 
as functions of n in the Matsubara frequency u>„ = 7rfcsT(2n4- 
1) for T/T c = 0.1, h /E F = 0.1, di = d 2 = 50k' 1 , and three 
values of the misorientation angle: a = (solid curve), ir/2 
(dotted curve) and n (dash-dotted curve). The curves are 
calculated in the clean limit, r — » oo, and for <j> = and 
6 = 0. 



FIG. 8: Spatial dependence of singlet and triplet pair am- 
plitudes f 3 and /to, normalized to the bulk singlet pair am- 
plitude /sj,, for parallel magnetizations, a = 0, T/T c = 0.1, 
ho/EY = 0.1, di = d,2 — SOfcp 1 and for two values of the 
impurity-scattering rate h/2rEF = 0.001 (thin solid curve) 
and 0.01 (thick solid curve). 



in Figs. 2] and [5] for thin and thick ferromagnetic lay- 
ers, respectively, for 8 = and for = 0. Note that 
f s (9 = tt) = / s (0 = 0), and f tl (6 = tt) = / tl (0 = 0), 
while / t o (6> = 7r) = -/to(0 = 0). For = the singlet 
amplitude / s is real, while the triplet amplitudes fto and 
fti are imaginary (Figs. [4] and [5]), and for — K the 
opposite is true (Fig. [5]). It can be seen that both spin 
singlet and triplet pair amplitudes decay in the same os- 
cillatory manner with distance from the FS interfaces. 
Similar behavior of these correlations, i.e., the absence of 
monotonically decaying long-ranged triplet pair ampli- 
tude, which is in a striking contrast with the dirty-limit 
case^ explains the monotonic dependence of I c on a far 
from the — tt transitions. For a = 0, amplitudes f s 
and fto oscillate with x around zero and f t \ = 0. For 
a = 7r/2 the amplitude fto oscillates as a function of 
x as in the previous case, while f s and f t \ oscillate in 
the same manner above or bellow zero. We find that fto 
and fti oscillate in phase and with the same character- 
istic length (proportional to 1 /ho) as f s , regardless of a, 
when < a < tt. Phase of f s oscillations is shifted by 
tt. Wc also find that triplet correlations penetrate into 
the superconductors and monotonically decay over £g, 
the same characteristic length, on which the singlet am- 
plitude saturates. For a — tt, when the magnetizations 
are antiparallel no triplet amplitudes exist and f s mono- 
tonically decays away from the FS interfaces as in SNS 
junctions. 
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FIG. 9: (Color online) Spatial dependence of singlet and 
triplet pair amplitudes / s , fto and fti, normalized to the bulk 
singlet pair amplitude f s t, for orthogonal magnetizations, 
a = tt/2, T/T c = 0.1, h /E F = 0.1, di = d 2 = 50^ and for 
two values of the impurity-scattering rate H/2tEf = 0.001 
(thin solid curve) and 0.01 (thick solid curve). 
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FIG. 10: Spatial dependence of singlet and triplet pair am- 
plitudes f s and /to, normalized to the bulk singlet pair am- 
plitude f s b, for parallel magnetizations, a = 0, T/T c — 0.1, 
ho/Ev = 0.1, d\ — d,2 = SOOfcp 1 and for two values of the 
impurity-scattering rate Ti/2tEf = 0.001 (thin solid curve) 
and 0.01 (thick solid curve). 



In order to study the effect of finite r, we have cal- 
culated {g a a) in the clean limit (t — > oo) as the largest 
contribution to the term h(g)/2T in Eq. (|2.1| . In Fig. [7] 
the typical dependence of real and imaginary parts of 
(fl 'a a) on the Matsubara frequencies is illustrated for thin 
ferromagnetic layers and three values of a = 0, 7r/2,7r. 
Takeing R « sign(o;„) for any a, and neglecting con- 
tribution of / which is nontrivial only in the vicinity of 
uj n = 0, we calculate g from the clean limit solutions, 
Eqs. (|2.16[) - (|2.23[) . replacing only k, k ± with k, k ± , given 
by Eqs. (j235)l (IQ5]) . 




FIG. 12: (Color online) The current-phase relation I(4>) for 
weakly disordered ferromagnets H/2tEf = 0.001, for T/T c = 
0.1, ho/E-p = 0.1, di — d-2 = SOOfcp 1 , and hve values of the 
misorientation angle: a — (thick solid curve), n/3 (dotted 
curve), 7r/2 (thin solid curve), 2n/3 (dashed curve) and n 
(dash-dotted curve). 




FIG. 11: (Color online) Spatial dependence of singlet and 
triplet pair amplitudes f s , fto and fti, normalized to the bulk 
singlet pair amplitude f s b, for orthogonal magnetizations, a = 
it/2, T/T c = 0.1, h /E F = 0.1, di = d 2 = 500/fc" 1 and for two 
values of the impurity-scattering rate H/2tEf = 0.001 (thin 
solid curve) and 0.01 (thick solid curve). 
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FIG. 13: (Color online) The current-phase relation I(<j>) for 
weakly disordered ferromagnets H/2tEf = 0.01, for T/T c = 
0.1, Hq/Ef = 0.1, di — Gfo = SOOfcp 1 , and hve values of the 
misorientation angle: a — (thick solid curve), n/3 (dotted 
curve), 7r/2 (thin solid curve), 2tt/3 (dashed curve) and tt 
(dash-dotted curve). 

The influence of moderate disorder in ferromagnets 
on the pair amplitudes is illustrated in Figs. ISlfTTl for 
the parallel (a = 0) and the orthogonal (a = tt/2) 
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magnetizations, for two values of the scattering rates, 
H/2tEp — 0.001 and 0.01. It can be seen that the only 
difference from the r — > oo case is the faster decay with 
the distance from FS interfaces. However, all amplitudes, 
singlet and triplet, decay with the same characteristic 
length. We can not see that the triplets are longer- 
ranged than the singlets. The influence of scattering on 
impurities in ferromagnets on the Josephson current is 
shown in Fig. [12l for thick ferromagnetic layers and for 
h/2rEp — 0.001. Comparison with Fig. [3] shows that 
current-phase relation is almost sinusoidal and the crit- 
ical current is three times smaller. For ten times larger 
scattering rate, H/2tEf = 0.01, numerical calculation 
gives exactly sinusoidal I{4>) and four order of magni- 
tudes smaller critical current, Fig. [T3l 

IV. CONCLUSION 

We have studied the Josephson effect in clean and 
moderately disordered SFFS junctions containing con- 
ventional (s-wave) superconductors, two mono-domain 
ferromagnetic layers with arbitrary angle a between in- 
plane magnetizations, and fully transparent interfaces. 
We have calculated the Josephson current and both sin- 
glet and odd in frequency triplet amplitudes in the clean 
limit and for moderate disorder in ferromagnets by solv- 
ing the Eilenberger quasiclassical equations analytically 
in the step-wise (non self-consistent) approximation for 
the pair potential. Using quasiclassical approach we lim- 
ited our considerations to the case of equal Fermi wave 
vectors, relatively weak ferromagnets and transparent in- 
terfaces. For the Josephson current we reproduced previ- 
ous results for the clean limit 19 and have found crossover 



to the dirty limit case for finite impurity-scattering rate 
which leads to sinusoidal current-phase relation and con- 
siderably smaller value of the critical current. In partic- 
ular, transitions between and 7r states are induced by 
varying the relative orientation of magnetizations, both 
in clean and diffusive junctionsj 19 ' 28 i 29 However, in the 
clean limit this is simply the effect of decreasing influ- 
ence of the exchange potential with increasing misorien- 
tation of magnetizations and far from — it transitions 
the critical Josephson current monotonically depends on 
the angle between magnetizations. This can be seen in 
the calculated spatial variation of the pair amplitudes. 
We find that for < a < it, both spin singlet and triplet 
pair amplitudes in F layers decay in the same oscillatory 
manner with distance from the FS interfaces. This decay 
gets faster as the impurity-scattering rate in ferromag- 
nets is increased, but characteristic length of oscillations 
is unchanged. The computed triplet amplitude f t \ has 
opposite signs in the two magnet regions, penetrates into 
the superconductors and monotonically decays over the 
same distance, which is the superconducting coherence 
length, as the singlet amplitude f s saturates to the bulk 
value. In contrast to the dirty limit case^ no substantial 
impact of odd in frequency spin-triplet superconducting 
correlations on the Josephson current has been found in 
the ballistic and moderately diffusive regimes. 
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